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ABSTRACT 

It has recently been demonstrated that coherent relative streaming velocities of order 30kms _1 
between dark matter and gas permeated the universe on scales below a few Mpc directly after re- 
combination. We here use a series of high-resolution moving-mesh calculations to show that these 
supersonic motions significantly influence the virialization of the gas in minihalos, and delay the for- 
mation of the first stars. As the gas streams into minihalos with bulk velocities around lkms -1 at 
z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities 
of the halos, increasing the typical virial mass required for efficient cooling by a factor of three, and 
delaying Population III star formation by Az ~ 4. Since the distribution of the magnitude of the 
streaming velocities is narrowly peaked around a non- negligible value, this effect is important in most 
regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of 
the absolute number of minihalos that can be expected to cool and form Population III stars by up to 
an order of magnitude. We further find that the streaming velocities increase the turbulent velocity 
dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass 
function of the first stars. 

Subject headings: cosmology: theory — early universe — hydrodynamics — methods: numerical — 
stars: formation 



1. INTRODUCTION 

iTseliakhovich fc Hiratal (|2010f ) have recently shown 
that the acoustic oscillations before recombination give 
rise to streaming velocities between the dark matter 
(DM) and gas. Whereas the pressure of the photon- 
baryon fluid suppressed the growth of baryonic perturba- 
tions prior to recombination, the DM aquired streaming 



velocities of order v str , 



SOkms -1 at z = 1000 for 



typical la fluctuations. During recombination, the sound 
speed of the gas dropped from roughly c/y/3to 2 x 10 _5 c 
as it transformed from a radiation-dominated plasma to 
a largely neutral gas. As a result, the relative streaming 
velocities between dark matter and gas transitioned from 
a subsonic into a supersonic regime, with typical Mach 
numbers M c± 5. After this period, the gas remained 
thermally coupled to the cosmic microwave background 
(CMB) through Compton scattering with residual free 
electrons, so that the sound speed decayed as oc a -1 / 2 
until z ~ 200. Disregarding sourcing, the streaming ve- 
locities therefore decayed at most as ex a -1 , so that the 
Mach number dropped to M > 2 once the gas became 
fully adiabatic, and remained constant thereafter. 

Streaming velocities of this magnitude may alter the 
abundance of DM halos of mass ~ 10 6 M and lead 
to an increased bias and clustering on scales between 
a few and 100 Mpc, roughly corresponding to the Silk 
damping length and the sound horizon at recombina- 
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tion (ITseliakhovich fc Hiratal 120101 ). iDalal et al.l ([201 Of ) 
argued that the 21 cm signal before and d uring reion- 
ization could also be modified substantially. iMaio et all 
(2011) exploited the coherence of streaming velocities 
below scales of about 10 Mpc to investigate their influ- 
ence on the formation of the smallest bound objects. 
They found that the abundance of DM halos with virial 
masses 10 4 < M v i r < 10 7 M was reduced by a few per- 
cent, while the gas fractions of these halos were af fected 
by a somewhat larger amount. iStacv et al.l (|201Q[ ) used 
smaller box sizes and concentrated on the collapse of gas 
in minihalos, finding almost no increase in the virial mass 
and only a mild delay of collapse. 

In the present study we also focus on minihalos, where 
the effect is likely to be most pronounced. Compared to 
previous investigations, we employ a higher resolution, 
which allows us to investigate the influence of streaming 
velocities in more detail, but limits the number of real- 
izations that we can investigate. Furthermore, we here 
employ the moving mesh code AREPO, which is based 
on a set of gri d points that are advected with the flow 
(Springel l2Q10h . This allows a higher efficiency and ac- 
curacy compared to SPH simulations, while maintaining 
the natural adaptivity of SPH. Compared to AMR sim- 
ulations, the Galiliean-invariance of AREPO is advanta- 
geous when dealing with supersonic bulk velocities. 

The structure of our work is as follows: In Section 2, we 
present the numerical setup of the simulations, followed 
by our main results (Section 3) and a brief discussion of 
their implications (Section 4). All distances quoted in 
this paper are in proper units, unless noted otherwise. 

2. SIMULATION SETUP 

We perform three statistically independent simulations 
in boxes with a side length of 500 kpc, denoted MH- 
1 to MH-3. For each realization of the cosmological 
density field, we first run a coarse DM simulation to 
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Fig. 1. — A comparison of three statistically independent minihalos with no streaming velocity (top panels), and with an initial streaming 
velocity of 3kms -1 applied at z = 99 from left to right (middle and bottom panels). We show the density-squared weighted gas temperature 
projected along the line of sight when the hydrogen density in the center has just exceeded nu = 10 9 cm -3 (top and bottom panels), and 
when the streaming case has evolved to the same redshift as the no-streaming case (middle panels). In the presence of streaming velocities, 
the effective Jeans mass of the gas is increased. The underlying DM halo therefore becomes more massive before the gas can cool, which 
delays the onset of collapse. We also find that virial shocks are more pronounced in the direction of the incoming streaming flow than in 
other directions. Non-linear effects of this sort may result in a higher velocity disperion of the gas (see also Figure 4). 
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Fig. 2. — The gas fractions and central densities of the miniha- 
los as a function of redshift. The different colors distinguish the 
individual realizations. The solid and dashed lines denote cases 
with no streaming velocity and a streaming velocity of 3kms — 1 at 
z = 99, respectively. The thick dashed line denotes the cosmologi- 
cal baryon fraction. The increased effective Jeans mass of the gas 
leads to a reduction of the gas fractions by about 50%, which in 
turn results in a lower central density and a substantial delay of 
the collapse of the gas. 



identify the Lagrangian region of the first minihalo that 
exceeds a virial mass of 5 x 1O 5 M . In these simu- 
lations, we employ 256 3 particles of mass 272 M and 
use a comoving gravitational softening length of 68 pc. 
The initial fluctuation power spectrum is that of a A 
cold dark matter (ACDM) cosmology at a starting red- 
shift of z = 99, with matter density Q m = 1 — £7 a = 
0.27, baryon density — 0.046, Hubble parameter 
h = iVlOO km s^Mpc -1 = 0.71 (where H is the 
present Hubble expansion rate), and a spectral index of 
n s = 0.96. Thes e parameters are base d on the five-year 
WMAP results (jKomatsu et al.ll2009D . The three real- 
izations use different linearly extrapolated present-day 
normalizations erg, which we set to as = 0.9, 1.1 and 
1.3, respectively. This is equivalent to altering the red- 
shift of collapse, and allows us to extend the scope of our 
simulations in light of the limited box sizes employed. 




Fig. 3. — The virial masses and collapse redshift s of all miniha- 
los for no streaming velocities (lower symbols), and for an initial 
streaming velocity of 3kms -1 at z = 99 (upper symbols). As in- 
dicated by the arrows, the virial mass required for efficient cooling 
is typically increased by a factor of ~ 3, which delays Pop III star 
formation by Az ~ 4. 



Once the location of a minihalo in each realization 
has been identified, we use t he refinement procedure de- 
scribed in lGreif et al.l (|2Qllh to increase the mass resolu- 
tion within a sufficiently large region around the target 
halo. The final masses of the DM particles and the mesh- 
generating points that represent the gas particles are 3.53 
and 0.72 M , respectively. The comoving gravitational 
softening length of the refined DM component is 17 pc. 
We displace the gas particles relative to the DM parti- 
cles by half the mean interparticle separation, implicitly 
assuming that the transfer functions of the gas and DM 
components are equivalen t. This can subs tantially affect 
the filtering mass scale (jNaoz et al.ll2010} ). but we chose 
to neglect this caveat since we here investigate the differ- 
ences between simulations with and without streaming 
velocities, but otherwise identical parameters. 

In simulations MH-l-REL to MH-3-REL, the gas par- 
ticles obtain an additional velocity of ^stream = 3kms _1 
in the x-direction, which approximately corresponds to 
a 1<j peak. The high resolution region is extended in 
the x-direction to ensure that only refined particles col- 
lapse into the minihalos. Simulations without an ad- 
ditional streaming velocity are denoted MH-l-NOREL 
to MH-3-NOREL. The chemical and thermal evolution 
of the gas is captured with a primordial chemistry net- 
work that tracks the abundances and the cooling from 
H, H+, H-, H+, H?, He, He+, He++ D, D+, HP, an d 
free electrons (| Glover & Jappsen 2007; Greif et al. 2011). 
We account for all relevant cooling processes, including 
H2 and HD rotational and vibrational line cooling, H2 
collision-induced emission (CIE) cooling, H2 collisional 
dissociation cooling, and heating due to three-body H2 
formation. For further details see IGreif et al.l (l20TTf ). 

3. RESULTS 
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Fig. 4. — The temperature and velocity dispersion of the gas as 
a function of distance from the center of the halo, shown at a point 
in time when the hydrogen density has exceeded nn = 10 9 cm -3 . 
The different colors distinguish the individual realizations, while 
the solid and dashed lines denote cases with no streaming velocity 
and a lcr streaming velocity of 3kms _1 at z = 99, respectively. 
The black dashed lines denote the approximate virial radius. In 
cases with a streaming velocity, the increased virial mass leads to 
a higher virial temperature. This difference is reduced at smaller 
radii (higher densities), since the gas cools very efficiently during 
runaway collapse. The bottom panel shows that the velocity dis- 
persion remains systematically higher in cases with a streaming 
velocity. The increased amount of turbulence might later affect 
the fragmentation of the gas and alter the mass function of the 
first stars. 



3.1. Delay of Runaway Collapse 

To first order, streaming velocities may be considered 
as an additional source of pressure, which increases the 
Jeans mass of the gas: 



Mj 



"eff 



2G 3 /V /2 ' 



(1) 



with v^ff = Cg + ^stream- According to this simple ar- 
gument, th e minimum h alo m ass required for collapse is 
increased. iStacv et al.l (|2010l ) showed that this picture 
qualitatively describes their simulation results. The ef- 
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Fig. 5. — The comoving number density of minihalos that are 
expected to cool and form Pop III stars for no streaming velocity 
(solid line), and for an initial streaming velocity of 3kms — 1 at 
z = 99 (dotted line). The factor of ~ 3 increase in minimum virial 
mass leads to a reduction of the number of star-forming minihalos 
by up to an order of magnitude. The influence of Pop III stars 
on observables such as the 21 cm background or the reionization of 
the universe might therefore be substantially reduced. 



feet is particularly important in minihalos, since a la 
streaming velocity of lkms -1 at z ~ 20 is compara- 
ble to the sound speed of virialized gas in halos with 
M v i r ~ 1 .5 x 10 5 M^ - a lower limit for Pop III star for - 
mation (jMachacek et al.ll2QQHlO ? Shea fc Nor manll2 0071 ). 

In Figure 1, we show the temperature of the gas 
once the central hydrogen density first exceeds nn = 
10 9 cm -3 . We define this as the collapse redshift, since at 
this point in time the central gas cloud has evolved well 
beyond the initial Jeans instability and the density in- 
creases exponentially. In cases with a streaming velocity, 
the minihalos are more massive at the point of collapse, 
implying a higher virial temperature. The motions ac- 
companying their virialization also appear to be more 
complex. For example, the virial shocks in the direction 
of the incoming streaming flows are more pronounced 
than on the opposite side. This effect is particularly evi- 
dent in simulation MH-1, where a filament was by chance 
aligned perpendicular to the streaming velocity. 

In Figure 2, we show the gas fractions and maximum 
densities of the three minihalos as a function of time. In 
cases with an imposed streaming velocity, the gas frac- 
tions are reduced by about 50%. This agrees w e ll wit h 
the analytical predictions of iTseliakhovich et al.l ((2010). 
At a given point in time, the maximum density of the 
halo is reduced, so that the halo does not cool until it 
has accreted enough mass to compensate for this effect. 
The resulting shifts in virial mass and redshift are shown 
in Figure 3. The virial masses increase from 3.0, 3.1, and 
1.8 x 1O 5 M to 8.0, 9.4, and 6.2 x 1O 5 M , respectively, 
corresponding to factors of 2.7, 3.0, and 3.4. The red- 
shifts of collapse are shifted from 19.6, 22.7, and 31.7 to 
15.6, 19.5, and 26.1, resulting in delays of Az = 4, 3.2, 
and 5.5. We note that the factor of c± 3 increase in virial 
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mass is particularly important since it non-linearly af- 
fects the number density of DM halos that can cool and 
form Pop III stars (see Section 3.3). 

3.2. Thermal Evolution and Turbulence 

The increased momentum and energy input by stream- 
ing velocities also affects the thermal and turbulent evo- 
lution of the gas. In Figure 4, we compare the temper- 
ature and velocity dispersion within all three minihalos. 
In cases with a streaming velocity, the gas temperature 
near the virial radius is systematically higher. This is not 
suprising, since the virial temperature scales with virial 
mass and redshift as 



10 3 K 
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10 6 M 
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(2) 



A factor of ^ 3 increase in the virial mass and a com- 
paratively small decrease in redshift therefore result in a 
factor of ~ 2 higher virial temperature. Once the gas be- 
gins to cool efficiently, this difference is quickly reduced, 
however, and we find no clear difference remaining below 
r ^ 10 pc. 

In the bottom panel of Figure 4, we also show the 
velocity dispersion of the gas in radial bins after sub- 
tracting the bulk velocity of the halo and a spherically 
averaged radial infall and rotation velocity. Unlike the 
temperature, the velocity dispersion remains systemati- 
cally higher by a factor of a few even at very high den- 
sities. Evidently, the gas becomes more turbulent in the 
presence of streaming velocities. This could affect the 
mass function of the first stars, since a higher ratio of 
the velocity dispersion to the sound speed usually results 
i n increased fragmentation at later stages in the collapse 
([Clark et al.l 1201 ll ). As a result, the typical mass of a 
Pop III star may be reduced. 

3.3. Abundance of Minihalos 

How common is the above delay considering the mini- 
halo population as a whole? Since each component of the 
large-scale streaming velocity is distributed according to 
a Gaussian, the magnitu de of the velocity vector obeys 
a Maxwell distribution (Tseli akhovich et al.ll2010l ): 



Pstream = ^0\d) ^ 47TV. 



t exp 



, (3) 



where aid is the dispersion of each individual velocity 
component. The fraction of the universe with streaming 
velocities greater than 1.5 km s -1 at z = 99, which we 
consider a lower limit for the above delay to be signif- 
icant, may be found by integrating the above function 
from a/2 = aid a/3/2 to infinity, which yields approxi- 
mately 0.86. This shows that our results may be consid- 
ered representative for most of the volume. 

The cosmological number density of minihalos hosting 
Pop III stars may then b e estimated using the Sheth- 
Tormen (She th et al.l l2Q01) mass function: 



^mhO) 



n st (M 1 z)dM , 



(4) 



where we set M m i n = 1.5 x 10 5 M for the case of no 
streaming velocities and M m i n = 5 x 10 5 M for the case 



of a universal la streaming velocity, representing the fac- 
tor of ~ 3 increase in minimum halo mass. The resulting 
number densities should be considered upper limits, since 
not every halo at the low-mass end forms a Pop III star. 
We set M max = 1O 8 M , but note that our results are 
not sensititive to this parameter, since massive halos are 
rare. As shown in Figure 5, the number of minihalos 
that cool and form stars is reduced by up to an order of 
magnitude in the presence of streaming velocities. Such 
a large effect implies that streaming velocities should be 
taken into account when the influence of the first stars 
on observables is investigated. 

4. DISCUSSION 

We have found that supersonic streaming velocities be- 
tween the DM and gas substantially delay the onset of 
gravitational collapse in minihalos. The virial mass re- 
quired for efficient cooling is increased by a factor of ~ 3, 
which results in an average delay of Pop III star forma- 
tion by Az = 4. Streaming velocities also enhance the 
build-up of turbulence during runaway collapse, which 
could affect the fragmentation of the gas and hence the 
mass function of the first stars. 

Our results agree well with the simulations presented 
in lStacv et al.l (|2010l ) concerning the delay of collapse. In 
addition, and in contrast to their result, we find a signifi- 
cant increase in the minimum halo mass required for the 
gas to cool. This might be due to the ~ 10 times higher 
resolution employed in the present study, ensuring that 
the DM softening length is comparable to the size of the 
Jeans-unstable gas cloud. If this criterion is not fulfilled, 
the collapse of the g as might be artificially delayed (e.g., 
lYoshida et al.l 120031 ). We have investigated this for the 
case MH-1. We increased the initial DM and gas particle 
masses by factors of 8 and 64 to 28.2 and 225 M for the 
DM and 5.8 and 46.4 M for the gas. In these simulations 
we also decreased the number of cells enforce d per Jeans 
length from 128 to 64 and 32 (for details, see lGreif et al.l 
l2011f ). When decreasing the resolution, the virial mass 
of the minihalo at the point of collapse is increased from 
3.0 to 3.3 and 4.7 x 10 5 M , and the redshift of collapse 
is shifted from 19.6 to 19.2 and 18.1. On the other hand, 
an increase of the initial resolution by a factor of eight 
results only in a marginal decrease of the virial mass to 
2.9 x 1O 5 M , and shifts the redshift of collapse to 19.8. 
This shows that insufficient resolution can artificially de- 
lay the colla pse of the halo, w hich might be one of the 
reasons why iStacv et al.l (|201Q[ ) find significantly higher 
halo masses. 

Due to the large shift in mass in the presence of 
streaming velocities, the number density of minihalos 
expected to form Pop III stars is reduced by up to an 
order of magnitude, and is further modulated with a 
spatial coherence length corresponding to many Mpc in 
the present universe. This could leave a distinct im- 
print on the 21 cm back groun d and the reionization of 
the universe (jFurlanettd 12006b iBarkana fc Loebl 1200 If ), 
as w ell as affect the obser vability of Pop III supernovae 
(e.g. lBromm fc Lo eb 2002). Ultimately, this modulation 
might be inherited by the distribution of more massive 
objects forming at later times, thereby biasing large-scale 
structure estimates of cosmological parameters such as 
those related to the properties of dark energy. 
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